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Abstract: The partition function on the three-sphere of ABJM theory contains non- 
perturbative corrections which correspond to membrane instantons in M-theory. These correc- 
tions can be studied in the Fermi gas approach to the partition function, and they are encoded in 
a system of integral equations of the TBA type. We study a semiclassical or WKB expansion of 
this TBA system in the ABJM coupling k, which corresponds to the strong coupling expansion 
of the type IIA string. This allows us to study membrane instanton corrections in M-theory at 
high order in the WKB expansion. Using these WKB results, we verify the conjectures for the 
form of the one-instanton correction at finite k proposed recently by Hatsuda, Moriyama and 
Okuyama (HMO), which are in turn based on a conjectural cancellation of divergences between 
worldsheet instantons and membrane instantons. The HMO cancellation mechanism is important 
since it shows in a precise, quantitative way, that the perturbative genus expansion is radically 
insufficient at strong coupling, and that non-perturbative membrane effects are essential to make 
sense of the theory. We propose analytic expressions in k for the full two-membrane instanton 
correction and for higher-order non-perturbative terms, which pass many consistency checks and 
provide further evidence for the HMO mechanism. 
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1. Introduction 

The partition functions on the three-sphere of N > 2 Chern-Simons-matter theories can be 
reduced by localization to matrix models [1, 2, 3] which have been much studied in the last 
years. Surprisingly, these matrix models contain an enormous amount of information. When 
studied in the large N limit, this information can be decoded in terms of M-theory AdS duals. 
For example, the leading large N free energy can be seen to reproduce the gravity action evaluated 
on-shell, as first found in [4] in the case of ABJM theory [5] (see [6] for a review and a list of 
relevant references). 

Although most of the work done on these models has focused on the leading order contribu- 
tion at large N, from the point of view of M-theory and quantum gravity the most interesting 
information is contained in the subleading corrections. For example, the subleading logarithmic 
correction in N corresponds to a one-loop correction in quantum supergravity, as it has been 
shown in [7]. There are two types of corrections in the large N expansion: the perturbative 
corrections in 1/N, and the non-perturbative or exponentially small corrections. The pertur- 
bative corrections can be computed in closed form in a large class of J\f = 3 theories [8] , and 
they can be resummed in terms of an Airy function, as first shown in [9] for ABJM theory. The 
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non-perturbative corrections can be in turn divided in two types: the ones due to worldsheet 
instanton corrections in the AdS dual, and the ones due to more general membrane instanton 
corrections in M-theory [10]. The appearance of worldsheet instanton corrections was anticipated 
in [11], in the case of ABJM theory. They appear naturally in the 't Hooft expansion of the ma- 
trix model, at strong 't Hooft coupling, and were determined in a systematic, recursive way in 
[4], in a weak string coupling expansion. 

The study of membrane instanton corrections is more challenging, since these are non- 
perturbative effects both in the string dual and in the large N matrix model. In [12] some 
information about these instantons (like their action) was obtained from the study of the large 
genus asymptotics of the 't Hooft expansion, but no concrete recipe was given to calculate them. 
In [8], a new method was introduced to study the ABJM matrix model and its close cousins, 
based on an equivalence with a quantum Fermi gas. In the Fermi gas approach, the Planck 
constant is naturally identified with the inverse string coupling, and the semiclassical limit of the 
gas corresponds to the strong string coupling limit. One of the main virtues of the Fermi gas ap- 
proach is that it makes possible to calculate membrane instanton effects systematically, at least in 
the WKB expansion. This opened the way for a quantitative determination of non-perturbative 
effects in the M-theory duals to Chern-Simons-matter theories. 

In a recent paper [13], Hatsuda, Moriyama and Okuyama (HMO) made various crucial 
observations on the structure of non-perturbative corrections in ABJM theory, where the inverse 
string coupling is essentially given by the Chern-Simons level k. First of all, they noticed that 
the worldsheet instanton contributions to the free energy can be resummed at finite k by using 
a Gopakumar-Vafa representation. The resulting expressions are divergent for integer k 1 . Since 
the free energy is finite for any value of k, these divergences have to disappear in the final answer, 
and [13] suggested that they cancel against similar divergences in the contributions of membrane 
instantons, in such a way that the total sum of all non-perturbative effects at integer k is finite. 
We will call this the HMO cancellation mechanism. This mechanism is beautiful and natural, 
and we believe is of deep conceptual importance for the understanding of M-theory. It shows, in a 
precise and quantitative way, that the genus expansion based on strings is essentially meaningless: 
in the non-perturbative completion of type II A string theory at finite, integer k through M-theory, 
only the combination of membrane instantons and worldsheet instantons makes sense. 

In some cases, the HMO mechanism gives a set of constraints for the membrane instanton 
corrections at finite k. Using these constraints, as well as the first few terms of the semiclassical 
expansion at small k obtained in [8], an expression for the one-instanton membrane correction 
was proposed in [13] which passes many consistency checks. It reproduces for example the low 
order, non-perturbative corrections to the free energy, for small integer values of k. 

The purpose of the present note is twofold. It was pointed out in [8] that the Hamiltonian 
problem appearing in the Fermi gas approach to ABJM theory can be studied with a pair of 
TBA equations first considered in [14, 15] and studied in detail in [16]. This observation was 
exploited in [17, 18, 13], where the TBA equations were solved for finite k but small chemical 
potential. This makes it possible to compute the partition function of ABJM theory for finite 
k and small N. On the other hand, the Fermi gas can be studied in the WKB approximation, 
at small k, and this provided many valuable insights into the problem [8]. Our first goal in 
this paper is to develop a semiclassical or WKB expansion directly in the TBA equations. This 
leads to an algorithm which calculates the grand potential of the ABJM model systematically, 
as a power series in k, and arbitrary chemical potential. This method is more powerful than the 

1 This had been already noticed in 2011 in unpublished work by the second author and Pavel Putrov. 
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original WKB expansion of the Fermi gas studied originally in [8], and one can easily push the 
computation to higher orders. 

Our second goal is to use this information to further explore the cancellation mechanism 
proposed in [13]. We verify that the expression for the membrane one-instanton correction pro- 
posed in [13] agrees with the WKB expansion to high order, and we propose analytic expressions 
for the full two-instanton correction and for some higher-order terms in the non-perturbative 
expansion. Our proposal for the two-instanton correction passes all the consistency tests, and 
in particular agrees with the results for the grand potential at k = 1,3 obtained in [13]. Our 
WKB results confirm then the HMO cancellation mechanism and are very helpful in obtaining 
conjectural expressions at finite k. 

The organization of the paper is as follows. In section 2 we review some general results on 
non-perturbative effects in ABJM theory. In section 3 we explain how the Fermi gas approach 
of [8] can be reformulated in terms of the TBA system considered by Al. Zamolodchikov in [16]. 
In section 4 we study in detail the WKB expansion of the TBA system and explain how it leads 
to the semiclassical expansion for the grand potential considered in [8]. In section 5 we give 
an application of the WKB expansion to the calculation of membrane instantons: we test the 
conjecture for one- membrane instanton corrections in [13], and we propose an exact formula for 
the two-membrane instanton contributions. Finally, we conclude with some open problems. In 
the first appendix we list some useful results for the integration of generalized hypergeometric 
functions, and in the second appendix we list WKB expansions for some low order contributions 
of membrane instantons. 



2. Perturbative and non-perturbative aspects of ABJM theory 



The quantity we will focus on in this paper is the partition function of ABJM theory on the 
three-sphere, Z(N, k), which is given by the matrix integral [1] 
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This matrix integral can be calculated in the 't Hooft expansion 
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by using techniques of matrix model theory and topological string theory [23, 4]. In particular, 
one can obtain explicit formulae for the genus g free energies appearing in the 1/N expansion 



F(X,g s ) = Y,g 2 r 2 F 9 (X), 

9=0 
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The genus g free energies F g (X) obtained in this way are exact interpolating functions, and they 
can be studied in various regimes of the 't Hooft coupling. When expanded at strong coupling, 
they have the structure 

F g (\) = FV(\) + F»P(\). (2.5) 
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The first term represents the perturbative contribution in a', while the second term is non- 
perturbative in a', 

F 9 np (A) ~ O (e~ 2nV ^j . (2.6) 

The type IIA dual of ABJM theory involves the space AdS4 x CP 3 [5], and this geometry sup- 
ports worldsheet instantons wrapping a CP 1 C CP 3 [11]. The non-perturbative piece (2.6) was 
interpreted in [4] as the contribution of these worldsheet instantons. 

Besides the non-perturbative effects in a' , one can use the connection between the large-order 
behavior of perturbation theory and instantons to deduce the structure of non-perturbative effects 
in the string coupling constant. In [12] a detailed analysis showed that these effects would have 
the form 

exp (-fcvrv 7 ^) (2.7) 

at large A. It was also proposed in [12] that the source of these effects are D2-branes wrapped 
around three-cycles in the target space. An appropriate, explicit family of generalized Lagrangian 
submanifolds with the topology of MP 3 C CP 3 was proposed as an explicit candidate for these 
cycles. We will refer to these non-perturbative effects as membrane instanton effects: they can 
be interpreted as M2-instantons in M-theory [10] , where the M2-brane wraps a three-cycle inside 
S 7 /Zfc which is the lift of the three-cycle in CP 3 . Notice that these membrane instanton effects 
are invisible in ordinary string perturbation theory. 

Further information on the membrane intantons can be obtained by using the Fermi gas 
approach introduced in [8]. In this approach, one first notices (see also [19]) that the matrix 
integral (2.1) can be written as 

This in turn can be interpreted as the canonical partition function of a one-dimensional Fermi 
gas with a non-trivial one-particle density matrix 

11 1 1 

P[X1,X2) ~ 2^ (2 cosh 1/2 (2 cosh f ) 1/2 2 cosh (^) " ( ^ 

The one-particle Hamiltonian H of this system is then defined as 

p = e~ H , (xi\p\x 2 ) = p(xi,x 2 ), (2.10) 

and the Planck constant is 

h = 2irk. (2.11) 

The semiclassical or WKB expansion is then around k = 0, and it corresponds to the strong 
string coupling expansion in the type IIA dual. 

The Fermi gas approach suggests to look instead to the grand partition function (see also 

[20]) 

oo 

Z(n,k) = l+Y,Z(N,k)z N , (2.12) 

N=l 

where 

z = e fM (2.13) 
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plays the role of the fugacity and fi is the chemical potential. The grand potential is then defined 
as 

J( f i,k) = logZ( f i,k). (2.14) 
The canonical partition function is recovered from the grand-canonical potential as 

At large N, this integral can be computed by applying the saddle-point method to 

Z(N, k) = J d/i exp [J(/z, k) - fiN] . (2.16) 

As shown in [8], the grand potential is the sum of a perturbative and a non-perturbative 
piece, 

J(n,k) = JV(n,k) + J n P{ii,k). (2.17) 
The perturbative piece is a cubic polynomial in \x: 

jP^k) = < =^-fi 3 + B(k) f i + A(k). (2.18) 

The coefficients C(k), B(k) where computed in [8] for ABJM theory (in fact, the analogs of these 
coefficients for other N = 3 theories can be also computed in closed form). The coefficient A(k) 
can be computed in a WKB expansion around k = [8], and the all-orders result was conjectured 
in [21]. In this paper we will be interested in the non-perturbative piece J np (n,k). In general 
one has the following result: 

^"(M, k) = J M2 (/x, k) + J ws Gu, *;) + ■■■. (2.19) 

Let us explain this structure in more detail. The first term J M2 ([x, k) has the following expansion 
for fx » 1, 

J M V, k) = M^V 2 + b e (k)n + c e (k)) e" 2 ^. (2.20) 

This type of contributions to the grand potential was first found in [8] and it was interpreted there 
as due to membrane instantons, i.e. to M2-branes wrapping a three-cycle A4 C S 7 /Z^ which is 
a lift of a three-cycle in CP 3 . The positive integer t is the winding number of the wrapping. We 
will refer to the £-th term in the infinite series (2.20) as the the contribution of the ^-membrane 
instanton. The coefficients ae(k), be(k) and ce(k) are non-trivial functions of k. Using the WKB 
expansion of the quantum Fermi gas, it was shown in [8] that they have a perturbative expansion 
around k = of the form 

1 oo 

a e (k) = -^ae, n k 2n . (2.21) 

n=0 

A similar expansion holds for b# (k) and cg(k), with coefficients b^ n , c^ n , respectively. The WKB 
expansion was obtained in [8] up to order n = 2. One goal of this paper will be to obtain a more 
efficient method to calculate the WKB expansion. 

The second term in (2.19), J ws (fi, k), is the contribution of worldsheet instantons wrapping 
CP 1 C CP 3 . It has the following expansion for ji » 1, 



J ws (tx,k) = d m (k)e-^. (2.22) 

m=l 
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In [13] a very useful formula was proposed for J (/j,, k): by using the fact [23, 4] that the ABJM 
matrix integral is dual to topological string theory on local P 1 xP 1 , one can write 




(2.23) 



In this formula, n 9 d is the weighted sum of the Gopakumar-Vafa invariants [22] of local P 1 x P 1 , 

n9 d = E <*■ ( 2 - 24 ) 

di+d,2=d 




Figure 1: From the M-theory point of view, there are two types of non-perturbative effects in ABJM 
theory: M2-branes can wrap a cycle A4 C § 7 /Zfe which descends to an MP 3 C CP 3 cycle in the type 
IIA target; or they can wrap a cycle W which descends to CP 1 C CP 3 . The most general M2-brane 
configuration wraps I times the cycle M. and m times the cycle W. 

Notice that, from the point of view of M-theory, both types of instantons are due to M2- 
branes wrapping three-cycles: membrane instantons correspond to a three-cycle Ad which de- 
scends to a three-cycle wrapped by D2-branes in the type IIA target (if the proposal of [12] is 
correct, this three-cycle is an MP 3 C CP 3 ); worldsheet instantons correspond to M2-branes wrap- 
ping the three-cycle W = S 3 /Zfc, which descends to CP 1 C CP 3 . Therefore, from the point of 
view of M-theory, there are two types of three-cycles, Ad and W, and the most general M2-brane 
configuration wraps £ times the cycle Ad and m times the cycle W, see Fig. 1. We should then 
expect that the most general contribution to J np (/x) is of the form 



tm>l 



2£ + 



4m \ 



(2.25) 



The contribution J M2 (/i, k) corresponds to m = 0, while the contribution J ws (//, k) corresponds 
to I = 0. The existence of such bound states, labeled by two integers (I, m), was first indicated in 
[13] in order to interpret their data on the grand potential. From the point of view of M-theory 
they are completely natural, since they correspond to the most general set of supersymmetric 
cycles in the geometry. 

Can we compute the contribution of these non-perturbative states? In the case of worldsheet 
instantons, the expressions (2.23) and (2.24) reduce the problem to the determination of the 
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Gopakumar-Vafa invariants of local P xP . This can be done in many different ways, and there 
is no difficulty in calculating J ws ( / u, k) to any desired order. The calculation of J M2 ([i,k) can 
be done order by order in k and at all possible orders in the membrane instanton number by 
using the WKB method of [8]. In the remaining of the paper, we will develop a more efficient 
method, based on the TBA equations of [16], to calculate this perturbative expansion around 
k = 0. However, there is no known procedure to determine the contribution of bound states 
appearing in (2.25), not even at small k. This is probably the most important open problem in 
this subject, as we will mention in the final section of the paper. 



3. The Fermi gas/TBA approach to ABJM theory 
3.1 Predholm determinants and TBA 

We will now summarize some results from [16, 24] which are needed in this paper. 
In [16], the following integral kernel is considered 

This defines a homogeneous Fredholm integral equation of the second kind, of the form 

/oo 
K(e,9')f(6') = \f(6), (3.2) 
-oo 

Let A a be the possible eigenvalues appearing in (3.2), and let us introduce the Fredholm deter- 
minant 

E(z) = H(l + z\ a ). (3.3) 

n 

We will regard S as a grand canonical partition function. The grand potential is then given by 

J(z) = log E(z) (3.4) 

and it has the expansion 



where Zg is given by the integral 

r°° JL P -2«(0i) aq. 

Z *= IT —o e ( 3 - 6 ) 
J-oofJi 2cosh^±i 2vr 

with the periodicity condition 

0/+i = 0i. (3.7) 
We now introduce the iterated integral 

-2u(0i) 2u(d e ) 



/oo — -2u(t)i) lu(pn) 

^ cosh^cos h ^...eo 5h ^ '--- d *'- ^ (3 ' 8) 

R (d) = e- 2u{e \ (3.9) 



-oo wu.i 2 >,vuii 2 2 
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Notice that 



/oo 
d9R e (0) = (4ir) e+1 Z e+1 (3.10) 
-oo 

and the generating series 

R (e\z) = ^(-^) e Ri(0) (3.ii) 



satisfies 



It was conjectured in [16] and proved in [24] that the function R(9\z) can be obtained by 
using TBA-like equations which first appeared in the context of two-dimensional N = 2 theories 
[14, 15]. We first define 

R + (e\z) = l(R(e\z) + R(6\-z)), 



2 

R-(0\z) = ±{R(0\z)-R(0\-z)). 



(3.13) 



Let us now consider the TBA system 

2u{9) - e{9) + J ^ — cosHe _ 0l) , 



= -z f 
R+(0\z) = e- 6 ^ 



(3.14) 



2tt cosh(6' - 6»') ' 

Then, one has that 



/oo j /n 
-oo 7T 



arctan r]{6') ( 3 - 15 ) 



cosh 2 (^-^)' 

This conjecture has been proved in [24] for general u{9). In general, the system (3.14) has to 
be solved numerically, although an exact solution exists for u(9) = e e in terms of Airy functions 
[25]. 

3.2 From the ABJM Fermi gas to TBA 

Let us now find the relation between the Fermi gas approach to ABJM theory and the TBA 
system considered above (see also [18, 17]). We start from the expression (2.8) for the partition 
function. We can now regard the density matrix as an integral kernel with the structure of (3.1). 
The Fredholm determinant of this kernel is nothing but the grand partition function of the Fermi 
gas. We can then use the results of [16, 24] to write a TBA-like equation determining its grand 
potential. Let us consider the quantities 

e+1 x/iM-ii-x e ~ 2v(x) f°° dx 1 ---dx e e -Mxi) . . . e -M*t) 



p (x,x)-(x\p \x)- 27Tk J ^ {27Tky 2cosh ^y.. 2cosh ^ (3-16) 



where 



v(x) = ^ log (2 cosh |) . (3.17) 
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If we now change variables 

Xi = k6i. (3.18) 

and compare with (3.8), we find 

P e+ \x,x) = —L^x) (3-19) 

where we have denoted 

Rt(x) = Re(9=^y (3.20) 

The function Re(0) is calculated with the TBA system (3.14) and the potential 

1 / k9\ 
u(9) = - log (2 cosh y J (3.21) 

which depends explicitly on k. The grand potential is given by 

' dxR(x\z). (3.22) 



(3.23) 



dz 4irk 

Notice that the function R{x\z) can be written as 

„ , . . Airk I 
R(x\z) = ( x — — 

Z \ e H-» + 1 

The quantity appearing in (3.23) is essentially the full quantum-corrected version of the Fermi 
momentum pf(%) of the Fermi gas. Indeed, (3.23) can be computed by using the Wigner map, 
which associates to any quantum operator O a function Ow(x,p) in phase space (see [8] for 
details) : 

1 f°° / 1 \ 2 
R(x\z) = — dp -j =—Pf(x). (3.24) 



kzJ-oo Ve^-^ + l/w nz 
There is an important property of the TBA system of [16] which is worth discussing in some 
detail. Notice that the functions e(6), r](8) make it possible to calculate both R{x\z) and R(x\ — z). 
The last quantity is given by 

1 



rw , \ 47rfc 
R{x\ - z) = 



(3.25) 



and it corresponds to the same one-particle Hamiltonian (2.10) but with Bose-Einstein statistics. 
If we now take into account the expression (3.3), we deduce that for Bose-Einstein statistics there 
is a physical singularity at 

z = A 1 > (3.26) 

where Ao is the largest eigenvalue of the non-negative Hilbert-Schmidt operator p(xi,X2)- This 
singularity corresponds of course to the onset of Bose-Einstein condensation in the gas, and as a 
consequence the functions R±(x\z) will have singularities in the x-plane for z > Aq 1 . But this is 
precisely the regime in which we are interested, since large TV corresponds to fx S> 1. Of course, 
the singularity at large positive z cancels once one adds up R + and R . 

The appearance of this Bose-Einstein condensate has prevented the direct study of the large 
N limit of ABJM theory with the TBA equations (3.14), even numerically: the standard iteration 
of the integral equations does not converge when z is large enough. In the papers [17, 18] they 
study in fact the small z regime of (3.14) at finite k in order to extract the small z expansion 
of the grand partition function and therefore the canonical partition functions Z(N, k) for small 
N . The large N limit is then extracted from these small N results by numerical extrapolation 
[18, 13]. We will now propose another approach to study the TBA equations. 
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4. The semiclassical TBA system 



4.1 General aspects 

As we showed in the last section, the TBA system is equivalent to the Fermi gas picture. In 
particular, after setting 9 = x/k, the function R(x\z) is the Fermi momentum of the gas. As 
shown in [8] , one can do a systematic computation of all quantities in the Fermi gas in a WKB 
expansion, i.e. in a perturbative expansion in k. This suggests studying the TBA equation in a 
scaling regime where k is small and one can construct a perturbative expansion around k = 0. 
To do this, we set 6 = x/k and write (3.14) as 



re 

U(x) = e(x)+ / 

J — c 

7]{x) = —Z J 



dx' logfl + ^Qr')) 

2irk cosh i^rr-) 

k ' (4.1) 



oo dx , e -e(x') 



2nk cosh (^^) ' 
where 

U{x) = 2u(e=^j. (4.2) 

The second equation in (3.13) becomes 

„ , , , „ . . , f°° dx' arctan n(x') , . , 

R-(x\z) = R+(x\z) / — „ . /v (4.3) 

V ' ' V ' ' J-oc *k cosh 2 (^) V ' 

In the case of ABJM theory, where the potential is given by (3.21), we have 

U(x) = log ^2 cosh |) (4.4) 

and it is independent of k. Notice that in [16, 24] one considers a general potential u(9), with 
no k parameter to start with, but we can study a one-parameter deformation of the problem by 
considering a potential u(k0), in such a way that U(x) in (4.2) is independent of k. The original 
problem is then obtained when k = 1. 

The advantage of the /c-dependent equations (4.1), (4.3) is that they admit a systematic, 
perturbative expansion around k = 0, where they can be solved algebraically. This can be seen 
from the fact that, when k — > 0, the kernel becomes a (5-function: 

lim — = -5(x), (4.5) 

fe^0 2vrfccosh(f) 2 v h y 1 

and the integral equations become algebraic equations. The most convenient form of the TBA 
system for the small k expansion involves finite difference equations, and it was already considered 
(for k = 1) in [16] and specially in [24]. In more complicated models, the system of TBA integral 
equations, when written in terms of difference equations, is usually called the functional Y system. 
This form of the TBA system can be easily obtained by Fourier transform. If we denote 

F[f(x);p] = f(p) = 4= I" dx/(x)e^, (4.6) 

V27T J_oo 

we obtain, from the first equation in (4.1), 

2cosh (U(p)-?(p)) = F [log(l + r ? 2 (x));p] . (4.7) 
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The first factor in the l.h.s. of (4.7) is the displacement operator 



ikiv d ifc-TT d 

e 2 dx + e 2 dx. (4.8) 

Doing the same thing in the second equation of (4.1) we find that the two integral equations are 
equivalent (provided some mild analyticity conditions are satisfied, see [24]) to the two difference 
equations, 



< ( x + — j + e \x - — j = U \x + — J + U [x - — ) - log(l + 
Trik\ ( irik 



(4.9) 



(4.10) 



'/ I >' + ,y j ■ '/ l ^ .y j ' 

Equivalently, in terms of R + (x), we have 

1 + V 2 (x) = R+\x + —j R+ lx - —j expW fx + - + U (x - — 

_ . . / vrifc \ / 7ri/c 
-zi?+(x) = r ? lx+ — I + r? I x - — 

The equation (4.3) giving i?_ reads now, 

gl (* + g) _ gz (g ~ f) = 2ifc , 4 n) 

i?+(x + ^) i?+(x-^) l + r/ 2 (x)' V ' ' 

When A; = 1, these equations have been written down in [24]. 

The TBA equations, in the form (4.10), (4.11), can be solved systematically in an expansion 
around k = 0. Let us denote 

r(x) = R + (x), i(*) = |pM. (4.12) 



We introduce the ansatz 



r(x) = ^2r n (x)k 2n , 

n=0 

oo 

V(x) = Y,Vn(x)k 2n , (4.13) 

n=0 

oo 

t{x) = Y,Ux)k 2n . 

n=0 

We then make a Taylor expansion in k of the displaced functions in (4.10), and we solve order 
by order in k for r n , rj n . Once rj is known, t(x) can be obtained from 

t(x) = -^-tan- 1 (r ? ), C = ^ (4.14) 
7r sm (, Z 
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The solution for these functions at leading order in the k expansion is immediate, 



2 , 
t = --sin (£), 

7T 



^0 



(4.15) 



where we have introduced the variable 



t = Z ^~ U - (4-16) 



It is easy to see that, at each order in k 2n , we obtain from (4.10) two linear equations for r n and 
rj n , which are solved in terms of derivatives of lower r n /, rj n i with n' < n and of the potential 
U (x). The procedure can be easily automatized in a computer code to obtain the functions r n , fj n 
to any desired order. 

It is now straightforward to write down a power series expansion for R-(x), R(x). If we 
denote the coefficient of k 2n in this expansion as R- :Tl (x), R n (x), we find 



(4.17) 



R-, n (x) = r n (x)t (x) + B n (x), 
R n {x) = r n (x) (1 + t (x)) + B n (x), 

where 

n-1 

Bn(0 = Yl r m (0tn-m(0- (4-18) 
m=0 

At leading order we find, for example, 

Ro = r (l + t ) = — ^i==arccos(e). (4.19) 

Notice that both ro(x) and ijo(x) have a branch cut, as functions of £, along [l,oo). However, 
Ro(0 is holomorphic on the half-plane Re(£) > — 1, as it can be easily seen from the representa- 
tion 

mo = ^;-^= log (e + Ve^i) ■ (4.20) 



TTZ 



This is the semiclassical manifestation of the phenomenon we pointed out at the end of the 
previous section: the functions 77, R+ have a branch cut for 

z > 2 exp (min U) , (4.21) 

due to the appearance of a Bose-Einstein condensate. However, this branch cut cancels in the 
function Rq, since this function corresponds to Fermi statistics and there is no possible source of 
singularities for positive z in that case. 

We can test the result for Rq with a semiclassical calculation of (3.24) in the Fermi gas. At 
leading order in the WKB expansion, we have 
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where 

H(x,p) = log ^2 cosh + U(x), (4.23) 

and 

^ = ^/I d Vco sl l( i) + l - (4 ' 24 » 
which indeed can be explicitly computed and it agrees with (4.19) 

4.2 Semiclassical TBA for ABJM theory 

So far, we have considered the semiclassical TBA equations for a general potential U(x). In the 
case of ABJM theory, we have 

exp {U (, + f ) + U ( X - f ) } = ^ - 4sin* (£) (4.25) 

and it is useful to express everything in terms of £. It is then straightforward to study the 
semiclassical TBA system at higher order. For example, at next-to-leading order we find the 
equations 



z 2 



2710711 = 4 £2 [ 4 



7T 2 



(K)) 2 - r rg) + 2r on 



vr 2 r 2 



(4.26) 



with the solution 



-zn = 2r?i - —Tfo 

_ (z 2 (2g 2 + 3) - soe) 

7,1(6 32, 2 (l-a 7/2 



(4.27) 



It is easy to see from the recursion that the functions 

(i - ef n+i/2 r n (i), (i - ef n+i/2 vn(o, (i - e? n - l/2 uo (4.28) 

are polynomials in £. Moreover, we have found that the explicit solution for R n (0 has no 
singularity at £ = 1 and it is holomorphic in the half-plane Re(£) > —1, as expected from general 
principles. This involves a non-trivial cancellation of poles at £ = 1 in the sum of terms (4.17) 
giving i? n (£), and it can be regarded as a check of the procedure. 

Although we have not studied in detail the structure of the functions r n (£), r? n (£), some 
patterns can be easily observed. As we will see in a moment, it is useful to use, instead of the 
variable £, the variable u, defined by 

« = y ■ (4-29) 
The dominant terms in r n , % for z — > oo and u fixed seem to have the general structure 



r n (u,z) 

Vn(u,z) 



(l _ zhj?\ 3n+ ^ (2n)!2 14 ™- 2 ' 

1 " ' (4.30) 



z 6n-l u 6n+l ^2n 



i z 2 u 2 

r i6 



\3n+| (2n)!2 14n ' 
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This has been guessed by looking at the first orders in the explicit solution, but we have not 
found an analytic proof. As we will see in the next section, however, the above structure for r n 
implies the conjecture for the form of the one- membrane function a\{k) conjectured in [13]. 

4.3 The grand potential 

Our ultimate goal is to compute the WKB expansion of J M2 (//, k), therefore we have to compute 
the grand potential. This follows from (3.22). We will denote 

8 T 

■h = i~- (4.31) 

oz 

It is convenient to split this quantity w.r.t. the parity of z, 

8 I 

^ = JT + J~, (4.32) 



where 

After changing variables from x to 

u = ^ = sech (|) , (4.34) 

as defined in (4.29), we find 



Uf = - I — 7 =^= R± (zu/A) . (4.35) 



-I 

7T Jo 



1 du 



u 



The WKB expansion of R± (x) leads to the WKB expansion of the grand potential considered in 
[8], 

J z {^k) = Y,JzMk 2n ~ 1 - (4.36) 

n>0 

The integrals appearing in the calculation of J ZtH can be evaluated in terms of generalized hy- 
pergeometric functions. For J+ n , the answer involves the functions 2 -Pi, while for J~ n we also 
find the functions 3^2- The result can then be expanded at large z, and from this expansion, 
together with (2.20), one reads the small k expansion of a^(fc), b((k) and c^{k), for any i. In the 
Appendix we collect some useful results on hypergeometric integration which are needed in these 
calculations, as well as the resulting small k expansion of these coefficients for t = 1, 2, 3. 

As a simple example of this procedure, let us calculate the leading order correction to the 
grand potential, J 2) o- We have 

= n L TT^/T^^m = ^ K fe) ' (4 ' 37) 

and 

J7o = \ fdusm-U—) - I_ = _ *F 2 f 1,1,1; -,-;-V (4.38) 
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where we used (A. 5) and (A. 8). The above results for J^ agree with the calculation in [8]. We 
have also checked that the results for J 2) i and J 2) 2 agree with the results in [8]. 

It is interesting to notice that, in order to read off the coefficients ae(k), be(k), it is enough 
to calculate J+: by using the structural result (2.20) and assuming that the branch cut of the 
log is along the positive real axis, we find that 

jM2,+ = £ [(27rilogz - vr 2 ) £a e (k) + m(£b e (k) - a e (k))] z' 2 ^ 1 . (4.39) 

It is easy to verify that the leading order term in the expansion of (4.39) at large z comes from 
integrating the leading term of r n (u, z) in (4.30). By using (A. 3) we find 



L 



1 du u 6n ~ 2 lT(3n- 1/2 r(l/2 / 1 1 z 2 \ /AJ ^ 

- r = i ' { w 1 2 Fi 3n- -,3n+ -,3n; — . (4.40 

o Vi^ 2 (, _ z^\ 3n+ 2 2 r(3n) X V 2' 2' 'ley v ' 

V 16 ) 

At large z we have the logarithmic behavior, 

(*, - i, 3„ + i, 3„; |!) = r( _ 1/2 ^_ 1/2) ^-' log(-, 2 ) + ■ ■ ■ (4.41) 

If we put everything together and we compare the result with (4.39), we find that the ansatz 
(4.30) leads to the following perturbative expansion for the coefficient ai(k): 

4 ^ (-1)" /7rk\ 2n 4 /irk\ , d . 

?*gWr(T) =-s cos (t)- (4 - 42) 

This is precisely the result conjectured in [13] for this coefficient. 

5. Predictions for membrane instantons 
5.1 The HMO cancellation mechanism 

The Gopakumar-Vafa representation (2.23) of J ws (p,k) shows that it has double poles at all 
integer values of k. Since the original matrix integral is not singular for any value of k, there must 
be some way of canceling these divergences. The proposal of HMO in [13] is that, in the total 
non-perturbative grand potential, and order by order in e _M , there should be no divergences. 
This means that, in the sum (2.25) over all bound states (£, m) which contribute to a given order 
in e~^, singularities must cancel. This is the HMO cancellation mechanism. 

In general, since the contribution of generic bound states is not known, it is difficult to verify 
this cancellation mechanism in detail. However, some non-trivial information can be obtained 
by looking at low orders in the expansion. The case studied in detail in [13] is the case of one- 
membrane instanton, i.e. of contributions to the grand potential which go like e _2/J . Let us 
assume that k is an even, positive integer. Then, the only contributions to the term of order 
e _2/J come from 

(£,m) = (l,0) and (£, m) = (0, fc/2), (5.1) 

i.e. they involve only worldsheet instantons and membrane instantons. In this case, the HMO 
cancellation mechanism can be studied in detail. The worldsheet instanton poles at k = 2m are 
of the form 



d m (k)e- Am ^ k = {-l) m - 1 



\n> | 4Qz + 1) | 2/i 2 + 2 M + 1 | ^ (m) 



7T 2 (k — 2m) 2 7r 2 (k — 2m) mix 2 



e- 2 ^+--- (5.2) 
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where can be calculated from the Gopakumar-Vafa invariants of local P 1 xP 1 . According to 
the HMO mechanism, the poles have to be cancelled by similar poles in the membrane instanton 
contribution. Using these constraints, as well as the first three terms of the WKB expansion 
calculated in [8], HMO were able to propose an ansatz for the values of the coefficients a\(k), 
b\{k), c\(k) appearing in (2.20), 



a i( fc ) = -^ cos 



irk 
~2 



bi{k) 



ci(k) = 



cos 



irk 

T 



CSC 



irk 

T 



(5.3) 



2 5/t k o 
-3fe + 12 + 2 CSC 



irk 

T 



H — cot 

7T 



irk\ 



cos 



7rk 
~2 



One can check that these expressions have the right singularity structure to cancel the divergences 
in (5.2): near k = 2m, one has 



bi(k) 



4(-l 



\in— 1 



n 2 (k - 2m) 



+ 0(k-2m), 



Am 



4 1 
2m) 3m 



2m 
~3~ 



(5.4) 



+ 0(k-2m). 



n 2 (k-2m) 2 ir 2 (k 

Once the membrane instanton and the worldsheet instanton contributions are added for even 
integer k = 2m, the poles cancel and we find a finite contribution 



(-1) 



m—l 



A/i 2 + 2/1 + 1 



+ + 



mn^ 



3m 



2m 

If 



-2/x 



(5.5) 



This result reproduces a numerical calculation of the terms of order e~ 2tl in J(/i, k) for k = 2, 4, 6, 
and confirms the ansatz (5.3). 

The conjecture (5.3) shows that there is a simple, hidden structure in the contribution of 
membrane instantons to the grand potential. Thanks to our semiclassical TBA expansion, we 
can now compute the series expansion of these coefficients to higher order in k, as listed in the 
Appendix. Our results confirm the expressions in (5.3). In the case of a±(k), as we showed in 
the last section, this follows from the ansatz (4.30). 

5.2 The second membrane instanton 

In order to provide further tests of the HMO cancellation mechanism with the current techniques, 
we have to look at exponentially small terms which only receive contributions from states of the 
form (£, 0) and of the form (0, m). As we saw in the previous section, following [13], this is what 
happens for the terms ~ e~ 2 ^ when k is even. Similarly, when k is odd, the exponentially small 
term e _4/i has contributions only from the states 



(£, m) = (2, 0) and (£, m) = (0, k), 



(5.6) 



and it involves the two- membrane instanton. This case was only partially analyzed in [13]. They 
noted that the contribution of the worldsheet instanton for k = m odd is singular, and it has the 
behavior 

2/i 2 + n + 1/4 



d m (k)e 



-Amfi/k 



m 



_47T 2 (fc 



+ 



m) 



2ji + \ 

2ir 2 {k — m) 



+ 



+ v 



(m) 



mir^ 



e~ 4 » + 



(5.7) 
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The can be calculated from the Gopakumar-Vafa invariants and they are given, for m = 1, 3, 
by 



,(!) = 



,(3) 



(5.8) 



37 

3' ~9~' 

The poles appearing in (5.7) were also conjectured in [13] to be cancelled by similar poles in the 
membrane instanton contribution. Using this cancellation mechanism, as well as our data for 
the WKB expansion, we propose the following exact expressions for the coefficients of the second 
membrane instanton: 



a 2 (k) 



1 

TT 2 k 



(8 + 10 cos (irk)), 



4 1 
b 2 (k) = o~f~ (! + cos ( 7rfc )) + 7^ csc (nk) (17 + 24cos(tt/c) + 9 cos (2irk)) , 

TT Z k Z7T 



4 5cos(7r/c) 



49 . 7\ cot(vrA;) „, ,, 

— cos(7rA?) H — - + 5fccsc^(7rfc) 

24 6 / tt 



(5.9) 



+ COs(7r/c) 



5cot(7rfc) 21, 9/ 

— - + — kcsc 2 (irk) 

47T 4 



It is easy to check that the above expressions reproduce the WKB expansions presented in the 
Appendix. They also cancel the poles in (5.7). Indeed, one finds that, near k = m odd, 



b 2 (k) = - 
c 2 (k) = - 



1 



ir 2 (k — m] 
m 

m) 



+ 0(k-m), 
1 



A7r 2 (k-m) 2 2ir 2 (k 
The finite piece, after adding both contributions, is 

"4/i 2 + n + 1/4 



+ 



2m 



(5.10) 



m 



3m 



+ 0(k - m). 



+ z/ m ) + 



mir^ 



3m 



2m 

IT 



-4n 



For m = 1 , 3 we find 



4fi 2 + /i + 1/4 



7T^ 



e- 4 ^ and 



4^ 2 +^ + l/4 20 
3^ + ¥ 



-4/i 



(5.11) 



(5.12) 



respectively, which correctly reproduce the terms of order e~ 4fl in the expression for J(/x, k) with 
A; = 1,3 presented in [13]. Notice that, although the coefficient c 2 (k) is very different from the 
coefficient ci(k), their finite parts near the relevant poles have the same form. 

5.3 Higher membrane instantons 

The HMO cancellation mechanism can be studied in detail for one-membrane instantons and 
two-membrane instantons. The contribution of higher membrane instantons mixes with the 
contribution of generic bound states and it is difficult to extract information about the former 
without a more detailed knowledge of the latter. In particular, we don't have enough information 
about the singularity structure of the membrane coefficients as a function of k, and so far the 
only available information is contained in the WKB expansions. 

In the case of the ae(k) coefficients, however, it is easy to fit these WKB expansions to a sum 
of trigonometric functions. We find that, for odd (even) instanton number, the a^{k) are given 
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by a sum of cosines whose argument is an odd (respectively, even) integer times irk/ 2. Their 
explicit expressions, up to instanton number i = 5, are 



If fnk\ 124 /3vrfc\ f5nk 

a 3 (k) = - ^88 cos j + — cos ^— j + 4 cos 

a 4 (k) = ^-(364 + 560cos(vrfc) + 245 cos(2tt/c) + 48 cos(3vrfc) + 8 cos(4vrA;)), 

it k 

If fnk\ 4 f3irk\ 9104 /5vrA;\ /7vrA;\ 

a 5 (k) = (6080 cos ( — J + 4100 cos ( — ) + — ^— cos ( — J + 536 cos (—J 

f9irk\ /lbrfc\ /I3vrfc 

+136 cos — — + 24 cos — — + 4 cos 



(5.13) 

Finding a natural expression which fits our data for the expansion of the b#(k), cg(k) with t > 3 
is more challenging, and further information (including more data points in the WKB expansion) 
is probably needed. 



6. Conclusions and open problems 

In this paper we have developed a semiclassical approach to the TBA equations of [16]. When 
applied to ABJM theory, this reproduces in a more efficient way the WKB approach developed in 
[8]. Our results confirm the conjecture made in [13] for the one-instanton contribution. We have 
also proposed analytic expressions at finite k for the membrane instantons at order two, which 
are in perfect accord with the HMO cancellation mechanism. In addition, we have obtained some 
conjectural results at higher instanton number for the coefficients a^(fc). 

It is obvious from the results in this paper and its predecessors that there is a new and rich 
story concerning non-perturbative corrections in the M-theory dual to ABJM theory. Corrections 
to the partition function coming from intrinsic M-theory objects can be now computed in detail, 
and sometimes we can even guess their exact expression for arbitrary k. Moreover, the HMO 
cancellation mechanism shows very clearly that the genus expansion of type IIA string theory, 
although it can be resummed with a Gopakumar-Vafa representation, is essentially incomplete 
at strong coupling: only when the contribution of membrane instantons is taken into account do 
we find a finite answer for integer values of k (precisely the values for which we believe that the 
theory is defined non-perturbatively) . 

On the other hand, it is fair to say that we are only in a preliminary exploratory period of all 
these non-perturbative phenomena. Although the ABJM matrix integral and the TBA system 
contain detailed information about these instanton effects, it is not obvious how to extract it. The 
techniques proposed in [18, 17, 13] work for finite k, but they are based on a small z expansion 
which corresponds to small N. In this paper we have proposed a WKB expansion which is only 
valid for small k, but provides analytic large N results. Clearly, both approaches are insufficient, 
and they should be combined: what we need is the large fi expansion of the grand potential at 
finite k. This will very likely require a clever analysis of the TBA system which circumvents the 
problem with the Bose-Einstein condensate singularity pointed out in section 3. 

A method leading to results at large N and finite k would also determine the functions 
fe, m (fJ,,k) appearing in (2.25) for general bound states (maybe order by order in e" M ). In fact, 
it seems to us that the most important open problem is finding an approach to calculate the 
contribution of bound states, about which nothing is known. Even a WKB approach to the 



- 18 - 



problem would be useful. It was suggested in [8] that worldsheet instantons appear as quantum- 
mechanical instantons in the Fermi gas, and bound states might then appear as exponentially 
small corrections in fj, to a general quantum-mechanical instanton amplitude. This approach 
might give a first handle on this problem. Another possibility is to look for trans-series solutions 
of the difference equations (4.1), similar to what was done in [27] for the difference equations 
arising in matrix models. 

In addition to these general issues, there are a multitude of more concrete questions that come 
to our mind when we look at the conjectural expressions for the membrane instantons. Their 
contribution involves trigonometric functions whose coefficients are rational numbers and, very 
often, integer numbers. Is there some sort of Gopakumar-Vafa formula for them? Is there some 
duality between membrane instantons and worldsheet instantons? Can we compute some of these 
membrane instanton corrections directly in M-theory, by adapting for example the framework of 
[26]? 

Finally, it was pointed out in [8, 28] that any non-perturbative result in the context of the 
Fermi gas of AB JM theory can be interpreted as a non-perturbative result for topological strings 
in local P 1 xP 1 . In particular, the grand potential studied in this paper can be interpreted as the 
topological string free energy at large radius, and the membrane instantons computed here lead 
to corrections to this free energy of order 0(e~ 1 ^ gs ). Is there an interpretation for these effects 
directly in topological string theory? Are there "topological" membranes which lead to this sort 
of effect? 

As Bertolt Brecht would put it, "so many stories, so many questions." 



Acknowledgements 

We would like to thank Albrecht Klemm for useful conversations. This work is supported in part 
by the Fonds National Suisse, subsidies 200020-126817 and 200020-137523. 



A. Hypergeometric integrals 

In the calculation of the grand potential one finds two types of integrals. For J+, we have 
integrals of the form 

I 1 __du u m 

Jo vr^ 1 ( 1 _z^y s+1 *' 

After the change of variable 

u 2 = t, z 2 /16 = Z (A.2) 

they can be written as 

\[ * (1 - 0-^(1 - - T ~^f^ (s. + \, ^; 1 + = ; z) . (A.3) 

For J~, we have to do integrals of the form 

r du_«"-^-W, (A . 4) 

JO Vl - U 2 ( 2 2 M 2\ 3s +2 

1 16 
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Using that 

1 sin~ 1 (x) 
x Vl - x 2 
we can write the above integral as 



2*1 ( 1,1, 



du 



4 Jo vT 



,u m 1 



2 2\ ~3s 
2TZT x 



16 



^11,1; 5; 5^ 
' ' ' 2' 16 



(A.5) 



(A.6) 



To calculate this integral, we need two ingredients. First of all, we use the identity 132 in p. 436 
of [30] , with n = 1 (we correct a minor misprint in the statement of the identity) : 



(l-x)-^ 2 Fi 1,1,t:;z 



(1) 



3s 



(!) 



3s 



^la ! H;=; I Vi|'ff( 1 - I )-'-' 



(A.7) 



After using this identiy, and the change of variables (A. 2), the integral (A.6) can be reduced to 
an integral of the form 



f d ^ + i-i(i_i)&s + i-%>+i-i p F q ( ai ,--- ,a p ;h,--- ,b q ;Zt) 
Jo 

_ T(a p+ i)T(b q+1 - a p+1 ) 



(A.8) 



r(6,+i) 



p+\F q +i{ai, ■ ■ ■ , a p , a p+ i; b\, ■ ■ ■ ,b q ,b q+ \;Z) 



with p = 2, q = 1, as well as to a sum of integrals of the form (A.l) with half-integer s. 

The asymptotic behavior as z > 1 of the hypergeometric functions involved in these expres- 
sions can be easily found by using for example the Barnes representation. For the function 2*1, 
it is given by (see for example [29], eq. (9.7.7)): 



r( 7 )(-z) 



n-l 



r(7 — a)F(a + n) 



^ (»- ^ - l)!(a)i(l + z y 



i=0 



+ 



r(7 — a — n)T(a 



i=0 



(a + n)i(l - 7 + a + n) 



i!(n + i)\ 



[ip(i + 1) + ip(n + i + 1) 



(A.9) 



— if) (a + n + z) — ip(j — a — n — i) + log(— z)] z 1 n . 



B. WKB expansions 

Here we list the WKB expansion of the membrane instanton coefficients ai(k), bg(k), ce(k) for 
I = 1,2,3. 

kai(k) — „ + + ncon/on + 



^ 2 2 96 11520 2580480 928972800 490497638400 (Q ^ 

1 1 K 7 LJ1 + o (k 18 ) 

357082280755200 342798989524992000 v ; ' 
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18 , 2 5ir 2 k 4 Tr 4 k 6 vr 6 fc 8 ir s k 10 ir 10 k 12 ir 12 k 14 
ka 2 (k) = „- + 5r — — + — — - — — + o ^ ooon - , f7nnnign + 



7T 2 12 72 4032 362880 47900160 8717829120 , fi ^ 

2092278988800 v ; ' 

, />n 400 „ ,2 371^ 2 fc 4 11597r 4 fe 6 16373tt 6 A: 8 74113vr 8 A; 10 62408051^ 10 £i 12 

kai(k) = I- 70 A; I I 

v ' 3vr 2 24 576 92160 6635520 122624409600 



43949569vr 12 A 14 5465453813vr 14 /c 16 

H h O (A 18 ) . 

2550587719680 12242821054464000 v ; 

, . . 4 5A 2 67vr 2 A 4 19tt 4 A 6 247vr 6 fc 8 897r 8 fc 10 1430857tt 10 A 12 
kb x {k) = — - — + — — — — + + 100 ^ ,, nn ^ + 



vr 2 6 1440 48384 38707200 1226244096 669529276416000 

1637vr 12 A 14 118522319vr 14 A 16 ^/,ir\ 

_i I h CM At) 

30607052636160 87413742328872960000 v ; ' 

kh(k )_ 33 77k2 HSvr 2 ^ 4 1199vr 4 A 6 251^ 6 A 8 607^ 8 fc 10 3118137r 10 A; 12 
2( ' ~ tt 2 6~~ + 72 15120 + 120960 ~ 119750400 + 130767436800 

56863vr 12 fc 14 4740503tt 14 A; 16 

j I h Cm A; l 

261534873600 213412456857600 v ; ' 

,,„, 2560 l2 4769^ 2 fc 4 258689^ 4 fc 6 159091tt 6 A 8 262188523^ 8 fc 10 

kbi(k) = 188k ^ I 

v ' 9vr 2 ' ! 90 30240 172800 3832012800 

160190711489vr 10 A 12 15646173899vr 12 A 14 14718859878607vr 14 A 16 . 18 , 

+ 41845579776000 133905855283200 + 1365839723888640000 + ^ ' 



(B.3) 
(B.4) 



(B.5) 



, ,4 2\ 1 , 2 . l2 / 13vr 2 vr 4 \ 4 / 55vr 4 vr 6 \ 
fcci (fc) = -2 - - + - (vr 2 - 1) k 2 + -— - — A 4 + — — + 



(B.6) 

vr 2 37 ' 12 v " ' V 360 5767'" ' V 96768 ' 691207 

+ |--^ ^Wf — 3659 ' 8 )„» (B.7) 

' 38707200 15482880 7 \ 5573836800 122624409607 v ; 

( 7139277r10 _ ^ \ fc i2 , ffc 14X 
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